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Vacuum energy remains the simplest model of dark energy which could drive the accelerated ex- 
pansion of the Universe without necessarily introducing any new degrees of freedom. Inhomogeneous 
vacuum energy is necessarily interacting in general relativity. Although the four-velocity of vacuum 
energy is undefined, an interacting vacuum has an energy transfer and the vacuum energy defines 
a particular foliation of spacetime with spatially homogeneous vacuum energy in cosmological so- 
lutions. It is possible to give a consistent description of vacuum dynamics and in particular the 
relativistic equations of motion for inhomogeneous perturbations given a covariant prescription for 
the vacuum energy, or equivalently the energy transfer four-vector, and we construct gauge- invariant 
vacuum perturbations. We show that any dark energy cosmology can be decomposed into an in- 
teracting vacuum+matter cosmology whose inhomogeneous perturbations obey simple first-order 
equations. 
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The apparent acceleration of the Universe today has led cosmologists to consider a great variety of different 
models for the dark energy which apparently dominates the expansion. Vacuum energy is the simplest model 
of dark energy. If spacetime retains a non-zero energy density even in the absence of any particles then this 
energy density would be undiluted by the cosmological expansion and could drive an accelerated expansion 
as the density of ordinary matter or radiation becomes sub-dominant. The simplest example of constant 
vacuum energy density, V, in Einstein gravity is equivalent to a cosmological constant, A = SttGnV, and the 
discrepancy between the value of the energy density required by current observations |1| and the typical energy 
scales predicted by particle physics is the long-standing cosmological constant problem [2| ■ 

The idea of decaying vacuum energy is a recurring concept in attempts to explain the present acceleration 
[3-[ll|. Unlike other models of dark energy, vacuum energy does not introduce any new dynamical degrees 
of freedom. It is the energy density of nothing. There seems to be little to be learnt from simply assuming 
an arbitrary time-dependent vacuum to obtain the desired cosmological solution. Ideally one should have 
a physical model from which one could not only derive the time-dependent solution but also study other 
physical effects, including the evolution of inhomogeneous perturbations which can be tested against other 
cosmological observations. For example, vacuum fluctuations of free fields can support an averaged density 
proportional to the fourth-power of the Hubble expansion, V oc H* [12]. Such a vacuum energy would not 
in itself support an accelerated expansioUjbut other forms such as V oc H have been proposed fl3| better 
able to match the observational data [1Jt3- It is not clear whether or how such time-dependent vacuum 
models can be compared with observations in an inhomogeneous universe. In this paper we will argue that it 
is possible to give a consistent description of vacuum dynamics, and in particular the relativistic equations of 
motion for inhomogeneous perturbations, given a covariant, physical prescription for the local vacuum energy 
or, equivalently, the vacuum energy transfer 4-vector, Q^ = —V^V. Thus one should be able to subject vacuum 
models to observational constraints even in the absence of a Lagrangian derivation or microphysical description. 

One class of dark energy models that do have a well-defined microphysical prescription are quintessence models 
where the vacuum energy is determined by the self-interaction potential energy of one (or more) scalar fields, 
V{ip) |19l42ll |. A canonical scalar field has both potential and kinetic energy, therefore the vacuum potential 
energy interacts with the kinetic energy of the scalar field [22|. The introduction of a scalar field leads to a 
new degree of freedom and one can recover the Klein-Gordon equations of motion for the field from energy 
conservation, allowing for the energy transfer between the scalar field and the vacuum. Quintessence has been 
generalised to consider scalar fields also interacting with matter, known as coupled quintessence |23l - l28J . 

In this paper we will consider a general form for the vacuum energy which may directly interact with existing 
fields or fiuids without invoking any other additional degrees of freedom. We show that any homogeneous dark 
energy cosmology can be decomposed into an interacting vacuum-j-matter cosmology. By lifting this spatially- 
homogeneous solution to a covariant interaction one can study inhomogeneous perturbations which obey coupled 
first-order equations of motions for the matter density and velocity. We present the equations of motion for 
linear perturbations in Einstein gravity and identify gauge-invariant variables for the vacuum perturbations. As 
an example we consider perturbations of a generalised Chaplygin gas cosmology, decomposed into an interacting 
vacuum-|- matter. 



II. VACUUM ENERGY-MOMENTUM 

We define a vacuum energy V to have an energy-momentum tensor proportional to the metric: 

Ti; = -Vgi; . (1) 

By comparison with the energy-momentum tensor of a perfect fluid 

Ti; = Pg'i + {p + P)u^'u, , (2) 

we identify the vacuum energy density and pressure with p — ~P — V . Since p + P = the four- velocity of the 
vacuum, w^, is undefined, as expected in the absence of any particle flow. 

Note that the energy density and four- velocity of a fluid can be identified with the eigenvalue and eigenvector 
of the energy-momentum tensor ^ : 

Tj^u" = -pw'' . (3) 

Because the vacuum energy- momentum tensor ([T]) is proportional to the metric tensor, any four- velocity, u'^, is 
an eigenvector 

t>'' = -yu'^ Vm^, (4) 

and all observers see the same vacuum energy density, V, i.e., the vacuum energy is boost invariant. 
However there may be an energy transfer associated with the vacuum energy. Let 

V^t,^ = Q, . (5) 

For the vacuum energy-momentum tensor ([T]), since the metric tensor is covariantly conserved, we have 

Q, = -V,V. (6) 

We see that if the vacuum energy is covariantly conserved, Qi, = 0, then the vacuum energy must be homoge- 
neous in spacetime, i.e., V^y = 0. A non-interacting vacuum energy is equivalent to a cosmological constant 
in Einstein gravity. Conversely, an interacting vacuum, Q^ ^ 0, is inhomogeneous in spacetime. 

Although the vacuum does not have a unique four- velocity, we can use the energy flow, Q,y, to define a 
preferred unit four- vector in an interacting vacuum 

"'' " |v,^v-F|i/2 ■ ('^) 

normalised such that u^u'^ = ±1 for a spacelike or timelike fiow. Note however that the u^ defines a potential 
flow, i.e, with vanishing vorticity: V^^m^j = 

Assuming that the total energy-momentum for ordinary matter plus the vacuum, T^^ + I^', is covariantly 
conserved (as is required by the Ricci identity in Einstein gravity) we therefore have 

V;.T^' - -Q. . (8) 

One might consider a simple case where the energy transfer is proportional to the matter four- velocity, Q^ — Qu^ 
or equivalently u^ = u^, but in fact this is a restrictive assumption since the energy flow ([7]) is then a potential 
flow meaning that it cannot too admit any vorticity. In this case the matter velocity itself would be described 
by a scalar field, u^ oc V^y. 

III. FRW COSMOLOGY 

A. Spatially homogeneous background 

The Einstein equations for a spatially homogeneous and isotropic Friedmann-Robertson- Walker (FRW) cos- 
mology, with scale factor a and Hubble rate H — d/a, reduce to the Friedmann equation 

and the continuity equations for vacuum and matter 

p + 3H{p + P) = -Q, (10) 

V ^ Q. (11) 

In this case the vacuum and matter are both spatially homogeneous on spatial hypersurfaces orthogonal to 
the matter four- velocity, u^. The vacuum energy is undiluted by the cosmological expansion, but does have 
a time-dependent density in the presence of a non-zero energy transfer, Q ^ 0, where Q^, = Qu^. Note that 
the energy fiow, u'^, and matter velocity, m^, necessarily coincide in FRW cosmology due to the assumption of 
isotropy. 



B. Linear perturbations 

For simplicity we will consider linear, scalar perturbations about a spatially flat {K = 0) FRW metric [22ll2' 
[31j . with the line element: 

ds'^ = -(1 + 2(j))df + 2ad,Bdtdx' + a^ [(1 - 21/-)% + 2d^d^E\ dx'dx^ . (12) 

The energy and pressure of matter is given hy p + Sp and P + SP, and the four-velocity of matter is given by 

uf'=[l-(l),a-^d'v] , Uf,^[-l-cj),d^e] . (13) 

where we define d^v = a{dx'^ /dt) and 9 — a{v + B). 

Once we allow for deviations from homogeneity in the matter and metric, we must also allow for inhomogeneity 
in an interacting vacuum. As remarked earlier, the vacuum has an energy density and pressure, but no unique 
velocity. In particular the momentum of vacuum vanishes in any frame, {p + P)9 = 0. On the other hand the 
energy flow u defined in Eq. ([7]), can be written in analogy with the fluid velocity (|13p as 

u^ = [1 - , a-^d'v] , u,, = [-1 - , d,e\ . (14) 

where from Eq. ([7]) we identify 9 = —SV/V. 

The energy continuity equations for matter and vacuum become 

V^ / • \ 

dp + 'SH{Sp + SP)-'S{p + P)'ilj + {p + P)—[9 + a'^E~aB) = -5Q ^ Q<i) , (15) 



5V = 5Q + Q(f). (16) 

while the momentum conservation becomes 

{p + P)9-ic^^H{p + P)9 + {p + P)(t> + 5P = -f + clQe, 

-sv = f + Qe. (17) 

where the adiabatic sound speed c^ = P/p and following [22, [2^, |31| we decompose 

Q^^[-Q{l + ^)-6Q,d,if + Q9)]. (18) 

Note that the vacuum momentum conservation equation becomes a constraint equation which requires that the 
vacuum pressure gradient, di{—6V) is balanced by the force di{f + Q9). This determines the equal and opposite 
force exerted by the vacuum on the matter: 

- f ^6V + V9. (19) 

i.e., the fluid element feels the gradient of the vacuum potential energy. We see explicitly that perturbations of 
a fluid coupled to the vacuum with P = —p has no additional degrees of freedom, in contrast to a dark energy 
fluid with Px i= —px- Using the vacuum energy and momentum conservation equations to eliminate 5Q and / 
we obtain 

5p+'iH{5p + 5P)-?,{p + P)ij + [p + P)—[9 + a^E-aB^ = -5V , (20) 

{p + P)e-'iclH{p + P)9 + {p + P)(t, + 5P = 5V+{l + cl)Ve. (21) 

The metric perturbations obey the Einstein constraint equations 
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U + Hcf)) J i^ + H (a^E - aP) = -4ttG {6p + SV) , (22) 

iP + H(t)^ -ArrGip + P)9 . (23) 



As expected, the vacuum contributes to the energy constraint (1221) but not the momentum constraint (I23p . 

In summary, we have seven first-order scalar perturbations, (/>, tp, E, B, Sp, 9 and SV. We can eliminate two 
variables by choice of temporal and spatial gauge, and we have the two Einstein constraint equations ([2^ and 
(|23p . This leaves three variables to be determined from the two first-order energy and momentum conservation 
equations ((20|) and ([2T|). Thus we also need a physical model to describe the vacuum interactions, SV, in order 
to obtain a closed system of equations. 



C. Gauge invariant perturbations 

It is well known that metric and matter perturbations can be gauge-dependent. Under a first-order gauge 
transformation, such as t — > t + 5t{t,x'^), the fluid density and pressure transform as dp ^ Sp — p6t and 
SP — >■ 6P — PSt. However a quantity such as the non-adiabatic pressure perturbation, SPnad = SP — {P/p)6p, 
is gauge-invariant. 

Similarly the vacuum perturbation transforms as SV — > SV — QSt and 6Q — > 5Q — QSt. Given 5P = —5p = 

—6V and P = —p = —Q it is clear that the intrinsic non-adiabatic pressure perturbation of the vacuum is zero, 
as it is for any barotropic fluid with P = P{p). 

The fluid 3-momentum transforms as -^ 9-\-5t and the energy flow transforms similarly as — > + (5t so that 
the relative velcity 9 — 9 is gauge invariant. The comoving-orthogonal density perturbation (on hypersurfaces 
orthogonal to the fluid four- velocity u^) is thus given by the gauge- invariant combination 

<5pcom ^ Sp + p9. (24) 

Therefore the comoving orthogonal vacuum density perturbation is 

Spcom ^6V + V9. (25) 

This is in general non-zero, but the vacuum perturbation on hypersurfaces orthogonal to the energy flow u'^ is 
identically zero: 

Apcom ^SV + V9^0, (26) 

given 9 = —6V/V. In fact this is a consequence of the energy flow being the gradient of the vacuum energy, 
as given in Eq. ^, and therefore the orthogonal hypersurfaces are uniform vacuum energy hypersurfaces by 
construction. Note that we can therefore write the comoving vacuum density perturbation (j25p as 

<5/5com ^V{9-9) . (27) 

Therefore, if the energy flow follows the fluid four- velocity, u^ — u^, then we have 9 = 9 and the vacuum is 
spatially homogeneous on comoving-orthogonal hypersurfaces, Spcom = 0. 

The gauge-invariant vacuum density perturbation on spatially-flat hypersurfaces is given by 

C^-^J-jSV, (28) 

and this will in general be non-zero. Another gauge invariant expression for the vacuum density perturbation 
is the vacuum density perturbation on uniform-fluid density hypersurfaces, which describes a relative density 
perturbation 

^ = 3(C-C)^-3i7(f-|). (29) 

Therefore, if the vacuum energy is a function of the local matter density, V = V{p), then the relative density 
perturbation must vanish and the vacuum is spatially homoegeneous on uniform-density hypersurfaces, S = 0. 
The total non-adiabatic pressure perturbation due to any intrinsic non-adiabatic pressure of the fluid and the 
relative entropy perturbation between the vacuum and the fluid is then 

<5F„ad = SP-cJp+ QH^^p + P) ^ ■ ^^°^ 

This vanishes for adiabatic fluid perturbations, SP = c^Sp, and adiabatic vacuum fluctuations, 5* = 0, or a 
non-interacting vacuum, Q = 0. 

IV. DARK ENERGY MODELS 

A. Scalar field dynamics 

As remarked earlier, a self-interacting scalar field, (/3, has a vacuum energy, V{ip), which may be inhomogeneous 
if it interacts with the kinetic energy of the scalar field. Consider a single scalar field with Lagrangian density 
£ = F{X) - V{ip) where X = -g^'^'V ^,LpSI ^Lpjl ^^^. In an FRW cosmology, a homogeneous scalar field 
with vanishing potential {V = 0) has kinetic energy density pp = 2XFx — F and pressure Pp = F fM*. In 
particular for F{X) oc AT" we have a barotropic fiuid with constant equation of state wp = Pp/pp = (2n — 1)^"'^ 



[3a |. which equals one for a canonical field, F ^ X. The energy flow from kinetic to potential energy, Eq.©, is 
given by Q^ - -{dV/d^)V^,^ [22]. 

The energy continuity equation ([5]) for the scalar field kinetic energy requires 

V^Tj^ = {F.xUip + FxxV^<^V^X) V,(^ = -Q, . (31) 

Hence we recover the generalised Klein-Gordon equation |21l] 

dV 

FxD^ + FxxV^V^''X- — ^0. (32) 

B. Decomposing dark energy models 

Numerous dark energy cosmologies have been proposed in terms of exotic dark energy solutions pdc(a) or 
equations of state Pde(Pde) and/or adiabatic sound speed c^. Some of these models have been proposed as 
unified dark energy models |3q which can play the role of both dark matter (P^o ~ 0) and vacuum energy 
(Pde ~ — Pdc) at different stages of cosmological evolution. 

Any dark energy fluid energy-momentum tensor ^ with density pdc can be described by pressureless matter, 
with density p„i and velocity u^j = u'^, interacting with the vacuum, V, such that pdc = Pm + V. The 
corresponding matter and vacuum densities are given by 

Pm = Pdc + Pdc, V ^ -Pdc ■ (33) 

while the energy flow is Q^ = V^Pdc- In an FRW cosmology this corresponds to Q = —Pdc- One might choose 
to decompose a dark energy model pdc{o-) into any two interacting barotropic fluids such that pde = Pi + P2, 
but this would double the degrees of freedom in the model unless one of these two "fluids" is the vacuum. 

Linder and Scherrer [33 showed that any dark energy cosmology could be represented by a constant vacuum 
energy, V, and a fluid with an suitable equation of state, Wf — {Pde + V)/{pde — V). In this case the vacuum 
energy was non-interacting but the fluid required a specific equation of state, whereas in our case we will consider 
pressureless matter, but the vacuum then requires a specific interaction.^ These are both examples of so-called 
dark degeneracy [sl, [3^ . 

One much studied dark energy model is the generalised Chaplygin gas, defined by the barotropic equation of 
state [sllii 



PgCg = -Ap-c"g. (34) 

This leads to a cosmological solution for the density 

/ , ^ 1/(1+") 

PgCg=(A + Pa-^(i+")) . (35) 

This has the simple limiting behaviour PgCg oc a~^ as a — ?► and PgCg — > ^1/(1+") as a ^ +cxo, therefore this has 
been proposed as a unified dark matter model. However such models are strongly constrained by observations 
since the barotropic equation of state defines a sound speed for matter perturbations which only reproduces the 
successful ACDM model when a ^ [ill, \^ . 

The decomposition (j33p into pressureless matter interacting with the vacuum has previously been considered 
for the generalised Chaplygin gas by Bento et al [43] . In this case we have the FRW solution 

V^A(A + Ba-^^'+"A , (36) 



and hence 

A={p,n + VrV. (37) 

The form of the FRW solution suggests a simple interaction 



1 In practice one could choose any equation of state for the matter and a specific interaction such that pm + V = pde- An example 
is a quintessence model, where a canonical scalar field has sound speed equal to one, which reproduces any given dark energy 
cosmology for a suitably chosen potential. 
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In the matter or vacuum dominated limits this reduces to an interaction of the form Q oc Hp,n or Q oc HV 
studied, for example, by Barrow and Clifton 44]. 

It is intriguing to note that the FRW model can be defined in terms of an interaction p8D with a single 
dimensionless parameter, a, whereas when defined in terms of an equation of state (j34p its definition requires 
both a and the dimensional constant A which determines the late-time cosmological constant. In the interaction 
model, A, and therefore the late-time cosmological constant, emerges as an integration constant dependent on 
initial conditions. 

However, to study perturbations in the decomposed model we must "lift" the explicitly time-dependent FRW 
solution to a covariant model for the interaction. We appear to have (at least) two choices. 

Firstly we could require that the energy flow follows the matter four- velocity, Q^ = Qu^, such that 9 = 0. 
In any interacting matter -f vacuum model where the energy flow, Q^, is along the matter four- velocity, u^, the 
perturbation equations (|20)) and (|2T|) in an arbitrary gauge reduce to the equations for the matter density and 
velocity perturbation 

Sp^ + 3H6p,n-iprni^ + Prn—[0 + a^E^aBj = ^SV, (39) 

e + cj} = 0. (40) 

Because 9 = 9, the comoving density perturbation for the vacuum (j27|) is zero and hence there is no force 
exerted by the vacuum on the matter particles in Ea. (j2ip . There is no pressure perturbation in the comoving 
gauge and hence the sound speed for the interacting matter -|-vacuum is zero, c^ = 0. 

These equations are simplest in a comoving synchronous gauge where 9 = and cf) = 0. Eq. (j40p is automat- 
ically satisfied and the Einstein momentum constraint (j23l) then requires that ■0 = 0, i.e., tp = C{x^). Since we 
have 9 = in the comoving gauge, we also have 6V = 0. Thus the matter density perturbation equation p9p 
and Einstein energy constraint (I22p are the same as in the absence of any vacuum energy 

V^ / • \ 

Sp^ + 3H5p^ + p,n^-{a^E-aBj =0, (41) 

^P + H(a^E-aB] = ^■KG5p,n, (42) 

except that the background evolution is that for an interacting vacuum-f matter cosmology. 

These perturbation equations are strikingly simple, apparently the same as if we had treated the vacuum as 
spatially homogeneous from the start, but in fact these equations do allow a gauge-invariant vacuum perturba- 
tion, C, = — C(x*), and hence a relative entropy perturbation, S. From Eqs. (1^ and (|57|) we can identify 

S= (^"' + ^)' M (43) 

ap™ [prn + (1 + a)V) A' ^ ' 

Hence from Eq. ([M| we can identify the relative entropy perturbation with a change in the local equation of 
state and non-adiabatic pressure perturbation ([5(71). JPnad = —V5A/A. 

On the other hand the requirement that Q^ = Qu^ does impose a restriction on the type of matter field that 
can be coupled to the vacuum. It must be irrotational with u^ ex V^y. And we do not have a clear covariant 
model for the interaction Q = —u^Q^. By requiring 5V = —Q9, we construct the vacuum perturbation, 5V , 
iteratively given the matter velocity, 0, and the background energy transfer, Q in Eq. psp . 

Alternatively one might instead define a covariant model for the energy flow, Q^j, independent of the matter 
velocity. In this case we need to identify the homogeneous interaction ([55|1 as a total derivative, Q = V. Given 
the homogeneous solution (j37p the simplest covariant choice appears to be 

The comoving vacuum perturbation, Spcom in Eq. (j27p . is in general non-zero in this model since 9^9. 

In an FRW cosmology the interaction (|44p reduces to Eq. ([55]) as required ii A = 0. However the dimensional 
constant A now appears again in the covariant defintion of the model, defining the interaction (j44p . and hence 
the late-time cosmological constant. If A is a universal constant, Vf^A = 0, then we have V^y ex V ^pm and the 
relative entropy perturbation (1291) must vanish, S = 0. In this case the interacting matter -|-vacuum, including 
inhomogeneous perturbations, is exactly equivalent to the barotropic fluid model for the Chaplygin gas (134p 
with adiabatic sound speed, c^ = P/ p- 

V. CONCLUSIONS 

In this paper we have considered vacuum energy as a source of spacetime curvature in Einstein gravity. We 
have shown that an inhomogeneous vacuum energy implies an interacting vacuum. The vacuum energy is a 



scalar potential and thus defines a particular spacetime foliation which, in a cosmological solution, will be 
spatially homogeneous. 

We have shown that in order to consistently consider inhomogeneous perturbations in the presence of a vacuum 
energy one must specify a covariant, physical model for the vacuum energy. Simply specifying a particular time- 
dependence for the vacuum energy in an FRW cosmology is unsatisfactory if the vacuum energy is introduced 
solely to produce a particular time-dependence of the cosmological expansion. It is more satisfactory to specify a 
physical model for the vacuum energy, or equivalently its interactions, from which the time-dependent solutions 
can then be deduced. Even if the time-dependent solutions are no different, one can go on to consider the 
behaviour of inhomogeneous perturbations which can be compared against further observational constraints. 

The familiar example of an interacting vacuum is the potential energy of a scalar field which is inhomogeneous 
if the potential energy is transferred to the kinetic energy of the field. The equations presented here simply 
offer a different perspective on the dynamics of a scalar field with vacuum potential energy but do not change 
the dynamics. On the other hand our equations enable us to consider vacuum energy interacting with other 
forms of matter, including pressureless matter or radiation. We can therefore consider vacuum energy models 
which do not introduce any degrees of freedom beyond those already present in the cosmology, e.g., unified dark 
energy models. 

As an example, we have shown how the generalised Chaplygin gas cosmology can be re-derived as a solution to 
an interacting matter -|- vacuum model. In particular we have shown that the interaction can be defined by a single 
dimensionless constant and the late-time constant vacuum energy (which appears as a dimensional parameter 
in the original Chaplygin gas model) can instead be derived as a constant of integration in the matter -I- vacuum 
model. We leave it to future work to re-assess the cosmological constant problem in this particular model, 
and to study the constraints placed on the decomposed model from the full range of cosmological observations, 
including inhomogeneous perturbations and the growth of structure. 
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